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It has been predicted theoretically that cyclotron harmonic waves should propagate in warm 
plasma confined by a magnetic field, and show cutoff and resonance behavior associated with harmonics 
of the cyclotron frequency and the upper hybrid resonance frequency. In this paper, a simplified analy- 
sis based on the quasi-static approximation is presented and used to establish the relevant dispersion 
relations. Computer solutions of these are presented showing the influence of the magnetic field and 
discharge parameters on the propagation. A number of recent experimental results can be explained 
in terms of this type of wave motion and serve to establish its vahdity. This work is discussed briefly. 
The paper discusses some of the implications of this mode in laboratory and ionospheric plasma physics. 

1. Introduction 

Wherever possible in theoretical plasma wave propagation studies the analysis is simplified 
by neglecting thermal motions of the charged particles. This allow^s the plasma to be treated 
as a medium possessing a frequency-dependent, effective permittivity, which may be a tensor 
quantity if a static magnetic field is present. This approach is familiar from the literature of 
ionospheric propagation and may easily be extended to include the effects of collisions [Ratcliffe, 
1959; Budden, 1961]. 

It is recognized that the possible excitation of waves whose propagation depends on the exist- 
ence of nonzero thermal velocities is ignored in the cold plasma approximation. To study these, 
an approach based on the Boltzmann equation is required. This may take two directions. The 
first possibility, the hydrodynamic approach, is to work with moments of the Boltzmann equation, 
and leads directly to a description of the plasma behavior in terms of macroscopic quantities such 
as charge density, current, and electron pressure. Analyses of this nature have led to the pre- 
diction of acoustic waves: At high frequencies these have the well-known plane- wave dispersion 
relation 
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where propagation is as exp i{a)t — k' r) , (Op is the electron plasma frequency, k is Boltzmann's 
constant, Te is the electron temperature, and m is the electronic mass. Although (1) has never 
been verified experimentally, there is very strong evidence from scattering experiments carried 
out in the laboratory with inhomogeneous cylindrical plasmas, and described accurately by the 
hydrodynamic theory [Crawford, 1963, 1964a; Nickel et al., 1963; Parker et al., 1964] that the 
expression is a valid approximation. We shall take up this point again shortly. 

At low frequencies, electrostatic sound waves are predicted with phase velocity independent 
of charge density and frequency. For a plane wave, the dispersion relation is, 

--H-w) ' (2) 
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where M is the ion mass. These waves have been observed in laboratory discharges [Crawford, 
1961; Little, 1962]. 

The second method of employing the Boltzmann equation is to study individual particle tra- 
jectories, and then to integrate over velocity space to obtain macroscopic quantities. It was 
first pointed out by Landau [1946] that this procedure shows (1) to be incorrect. An additional 
damping term is required, due to trapped particles traveling with the wave [Bohm and Gross, 1949], 
which cannot be predicted by hydrodynamic theory. This so-called "Landau damping" becomes 
very strong at wavelengths comparable to a Debye length, and has been studied quantitatively 
recently in the laboratory [Malmberg and Wharton, 1964]. 

Since the single particle theory can be demonstrated to contribute damping terms to the 
acoustic modes predicted by hydrodynamic theory, a question now arises: does it also predict the 
existence of new types of propagating waves? The answer is. that it does. Cyclotron harmonic 
waves can propagate in warm plasma immersed in a magnetic field. The objects of this paper 
are to demonstrate physically how these waves can occur; to give computed dispersion curves 
for them, and to indicate the nature of the very recent experimental evidence confirming their 
existence. 

2. Theory 

2.L Single Particle Trajectory 

First consider the dynamics of a single electron moving under the influence of a static magnetic 
field and a transverse rf electric field, Eix, given by Ex exp i{cot — kxX — kzz) . We have, for the 
unperturbed trajectory of figure 1, 



Vr=V^ 



R 



Vz = V, 



(3) 



(f)' = ct)-\-a)c{t'-t). 



( TIME . t' ) 




Figure 1. Single particle trajectory (the unperturbed 
motion is along a helical path, the projection of which 
in the x-y plane is a circle). 
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where coc is the electron cyclotron frequency. Some simple integration and trigonometric manipu- 
lation, leads to 



Z' =Z-\-Vz{t' —t) 

:;c'=jc-f ^— j[sin {codt' -t) -^(f)} -s'ln (/)] 



(4) 
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For the jc-component of the perturbed trajectory, x[ , we have the equation of motion, 

x[-\-oy^[ = 'r)Ex exp i{o)t' — kxx' — k^z') , (5) 

where r] is the electronic charge-to-mass ratio. 

For small perturbations, the argument of the exponential may be replaced by its unperturbed 
value to obtain x[ correct to first order. This yields, 



x[ -hco^J =A exp i 



{w-k,v,){t'-t)-(^) sin {a>cU'-0 + (/>}j 



(6) 



A = r\Ex exp i 
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M'-c^y-*)-^' 



where we note that (vr/coc) is the Larmor radius of the particle motion. 

A well-known Bessel function identity may now be introduced [Bateman Manuscript Project, 
1953] 

exp -il^^y in {(Oc(t'-t)^(t>}] ^ 2 -^n (^) exp (-m{a>e(^'-0 + </)}). (7) 

Inserting this in (6) and integrating yields, 

x[=A exp [i{a) - k,v,){t' - 1)] V ^ "^^ , . j 7^. . (8) 

^f^^ [(ji)^ — i(x)-kzVz — n(x)cr\ 

The rf perturbation becomes very large when {oj — kzVz) approaches (/i±l)a)c, i.e., strong 
interactions occur between the rf electric field and the electron when the Doppler-shifted frequency 
approaches a cyclotron harmonic frequency. These cyclotron harmonic interactions result from 
two factors: the d-c rotational energy of the electron, and the fact that the rf electric field has a 
transverse variation (A:^- 7^ 0). This gives rise to space harmonics of the electric field over the 
direction of the electron motion, and hence strong interactions. As kx or vr decrease, the higher- 
order interactions will become progressively weaker. When {kxVRloc) is zero only the fundamental 
cyclotron resonance predicted by cold plasma theory remains. 



2.2. Multiparticle Treatment 

The next step is to develop a multiparticle treatment in which the space-charge forces due to 
the presence of many electrons couple their motions and lead to a dispersion relation. The first 
attempt to do this appears to have been due to Malmfors [1950], who studied the stabihty of crossed- 
field microwave devices. His work contained an algebraic error that was corrected by Gross [1951], 
whose work was in turn taken up and extended by Sen [1952], using the Laplace transform approach 
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initiated by Landau [1946]. The subject was discussed again in a series of papers, the most 
complete of these being that due to Bernstein [1958] which corrected several errors in previous 
analyses. 

In the averaging over the individual particle motions, two approaches are possible. A tensor 
expression for the rf current can be derived and substituted in Maxwell's equation [Bernstein, 
1958; Stix, 1962] or an approximation can be made at the outset which allows us to work with 
scalar quantities and simpHfies considerably the tedious integrations involved. This is the quasi- 
static, or "slow-wave," approximation, and consists in replacing the rf electric field E, by the 
gradient of a potential, (— V$). The approximation amounts to considering only waves for which 

\kr\'+\h\'>^\€ij\, (9) 

where eij are elements of the tensor permittivity of the plasma [Stix, 1962]. We assume, without 
loss of generality, variations in the x and z directions only, and c is the speed of light. It is also 
implied that E and k are parallel. Poisson's equation, for the rf space charge density is then, 

-iik.E, + k,E,) = {kl+kl)^ = ^- (10) 

To find a dispersion relation, it is now necessary to determine pi in terms of the rf electric 
fields in the system. The basic steps in the analysis will now be indicated without specifying the 
unperturbed electron velocity distribution. 

2.3. Determination of the RF Charge Density 

Liouville's theorem concerning particle number density, f{x, y, z, Vx, Vy, Vz) in phase space 
states that the number density of particles in a volume element is conserved as this volume moves 
in phase space. Hence, if /o and/i are the d-c and rf components of/, we have 

dfo.dfi 

rfT+dT^^- (11) 

Written out in full, this is 

'|^ + r-V/o + 7,irx5„.^] + [^+7r-V/:+T,(£+txfio+txfi,)-f|=0. (12) 

To obtain /i correct to first order, we must solve 



(f) 



\ 

where the time-differentiation is along the unperturbed electron trajectory, and the rf magnetic 
field term can be neglected in a quasi-static analysis. Integration of (13) gives 



fdt)=-vj'E'-^dt'+Mt,). (14) 

Next, assume that/i is of the form F exp i{o)t — kxX — kzz), and that E' is given by E exp i{(t)t' 
-kxx' — kzz'). Substitution from (3) and (4) in (14) yields in cylindrical coordinates, 
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F (v,,VR, <)>) = - 71 j' expi (a, - k,v,){t' - 1) - (^^y-m \^^{t'-t)j cos ^^ it' -t) + (t> 



E,^ + E,.^Ciys{coAt'-t) + <f>} 

OVz OVr 



dt', (15) 



A lower limit of ^o =~^ has been taken, assuming that there is no initial bunching of the electrons, 
i.e.,/i(^o)^ as Iq-^ — ^. This is equivalent to assuming that a small loss component is present. 
The limit is also justifiable as the only one yielding periodic solutions in 4> [Stix, 1962]. 

To obtain pi, F must be integrated over f^, vr, and (/). The integration over </) is between limits 
of and 2tt and may be accomplished by means of the substitutions of r for {t' — t) and i// for 
[(/) + (a>cT/2)], and use of the Bessel function identity [Bateman Manuscript Project, 1953] 



( 2kxVR 

Jn I SI 



(Oc 



^""2-j^2^Jo 



exp 



. /2kjcVR . oJcA . 



cos n^ dijj. 



(16) 



The resulting function, G, is given by 



ro r 

G{vz,vr) = — 2777) exp i{a) — kzVz)T 

J —oc 



r dfo J 2kxVR . COrT 

Ez-T—Jo sm — 



.^ a/o /0J,.t\ , l 2kxVR . OicT 



dr. (17) 



To complete the integration, a further Bessel function identity and its differential form are 
required [Bateman Manuscript Project, 1953] 
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(18) 



where Ao= 1; A„ = 2, n 9^ 0. Using these and again assuming that a slight loss component is pres- 
ent, we obtain 



G{vz,VR) = — 277r]<P 



\n{o)-kzVz)Jli-^^] 



k,^y 



\ Vr J dVR ^^ 



2 72 kxVR 



nV\ 



(w — A:jt;j)2 — (nwc)^ 



(19) 



An equivalent expression to this has been derived by Harris [1959, 1961]. Its value lies in the fact 
that when the unperturbed velocity distribution,/©, is specified, the rf charge density is immediately 
obtainable from 



Pi 



r +00 r 00 
= e I dvz I 

J -oc Jo 



vrG{vz, VR)dvR • 



(20) 



If several species of particles are present their space-charge contributions are additive. 
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2.4. Dispersion Relations for Oblique Propagation (kx9^0, kgy^O) 
It is instructive to consider first a uniform electron beam having a velocity distribution of the 



form 



/o = 



no 



2ttv, 



RO 



8{Vr — Vro) ^iVz — Vzo)' 



(21) 



The 8-functions imply that the beam has a translational velocity, v^o, and that all electrons have 
the same speed of gyration, vro, about the field lines. Substitution of (21) in (19) and (20), and 
combination of the result with (10) leads to the dispersion relation for the cyclotron harmonic 
waves on the beam. 



= kl 



1-^^S 



^'U''^^^) 



+ kl\l-wl'^AJli(ji) 

L n=0 



{(o — kzVzoy-^ {nojcy 



[{(o — kzVzo)^— {ncocy]^ 



}]• 



(22) 



where cop is the electron plasma frequency corresponding to the beam density, and jjl has been 
written for (kjcVRolcoc). It will be noted that this is the parameter previously encountered in the 
single-particle model of (8). As jjl^O we retrieve the Doppler- shifted cold-plasma dispersion 
relation 



= kl 



(co — kzVzo)^ — (o^ 



+ ^! 1- 



(co — kzVzof. 



(23) 



Next, consider an isotropic Maxwellian electron velocity distribution such that 



^'' = "»(2^)"""P 



-fn(r| + t;|) 

2KTe 



(24) 



Substitution of (24) in (19) and (20) and combination with (10) will again lead to a dispersion 
relation. Alternatively, (22) may be regarded as the basic dispersion expression which is then 
written with (24) as 
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dvz 



(25) 



where \ has been written for {k%KTelmo)^) or (kxR)^; and hence, R is the radius of gyration of the 
average thermal electron. 

The integration with respect to /x is accomplished using the identity [Bateman Manuscript 
Project, 1953] 



ri?(M)exp (-^) iidfJi = k exp (-X)WM. 
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(26) 



After some manipulation we arrive at the expression 



k%+ki=<oi(2^r 



^ exp ( ^^ j dvz 
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dvz 



(27) 



The integrand of (25) regarded as a function of vr (or jx) has no poles. Consequently, no Landau 
damping terms appear in (27) due to the electron motions transverse to the magnetic field. The 
physical reason for their absence is that all electrons retrace their paths periodically, at the cyclo- 
tron frequency. Hence, there is no phase-mixing. However, (27) indicates that there can be 
Landau damping due to the velocity distribution in the z-direction, except for the special case of 
purely perpendicular propagation {kz — O). Oblique propagation has been dealt with by Stix [1962] 
and we shall not discuss it further other than to remark that the damping will only by negligible 
if the condition {(t) — n(x)c)>kz(2KTelmy^'^ is satisfied. This implies that for modes corresponding 
to n> I propagation must be almost perpendicular to the magnetic field. By far the most inter- 
esting case, then, is that of perpendicular propagation which we shall now discuss in more detail. 



2.5. Dispersion Relations for Perpendicular Propagation (A^^O, k^ 
For the two velocity distributions studied, (22) and (27) give 
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It may be remarked here that the series representation is not the most satisfactory for computa- 
tion. An extremely useful integral expression is derived in an appendix to the paper. 

The behavior of the two expressions of (28) differ considerably. The first is shown plotted 
in figures 2a, b, c, and d for different values of (co^/co^). Each passband demonstrates a resonance 
(kx^^) at the appropriate cyclotron harmonic frequency. Cutoffs (A:j--^0) occur at (/i-f l)ajc, 
and at the point X corresponding to the upper hybrid frequency (co^ + co^)^/^ familiar from cold- 
plasma theory. The oscillatory nature of the dispersion relation in each passband stems from the 
nature of ^/^^t[y^(^)] which shows all of the zeros of Jnifi) and J'JilJL) as indicated in the figures. 
At low values of (cOp/co^) individual modes do not couple. P"or larger values, such coupling can 
occur and gaps in the propagation appear. In these gaps the solution for o) obtained for real 
values of /x is complex, and indicates the existence of growing waves. This implies that the sys- 
tem is unstable [Crawford and Tataronis, 1965]. We shall return later to this extremely important 
point. 

Figure 3 shows dispersion curves for the Maxwellian distribution. The same resonances 
and cutoffs are indicated as in figures 2a, b, c, and d but no coupling of the individual modes occurs 
and the system is absolutely stable. 

So far, we have not discussed the effect of collisions on the otherwise undamped perpendicular 
propagation of cyclotron harmonic waves. It has been shown elsewhere [Crawford, Kino, and 
Weiss, 1964b] that electron/neutral collisions occurring with frequency v can be taken into account 

in (28) by replacing co by (co — i^') and oj^ by co^ I 1 — i — J. To obtain a simple order of magnitude 

estimate for the attenuation, it may be assumed that kx is some function, K, of (o) — iv). Evaluating 
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Figure 2. Dispersion characteristics for perpendicular 
propagation of cyclotron harmonic waves — monoener- 
getic transverse velocity distribution. 
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Figure 3. Dispersion characteristics for perpendicular 
propagation of cyclotron harmonic waves — Maxwellian 
transverse velocity distribution. 




(k^R) 



the real and imaginary parts of /cj-, yields 



kxr I ikj 
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(29) 



The group velocity, Vg, may be introduced to write k^i in the form (— v/vg) which represents attenua- 
tion/unit length. It is a feature of the computed curves of figures 2a, b, c, d, and 3 that Vg generally 
lies below the electron thermal velocity and decreases rapidly with increasing harmonic number. 
This means that the collision damping of perpendicular propagation may be strong if v is not small, 
and will always limit the number of harmonics observable. 

3. Experiments 

3.1. General Comments 

In this section we shall describe some of the experiments which confirm the existence of 
cyclotron harmonic wave?. These are dealt with in detail elsewhere [Crawford, 1964b], and we 
shall only indicate briefly the nature of the phenomena involved. It should be emphasized that 
although the observations date back to 1959 [Wharton, 1960], their explanation was not then 
clear, and it was in 1963 [Canobbio and Croci, 1964; Tanaka, Mitani, and Kubo, 1963] that the 
suggestion was first advanced that these waves might be responsible for certain cyclotron harmonic 
noise phenomena in plasmas. Since these were the subject of the earliest studies we shall begin 
our brief review with some comments on them. 



3.2. Studies of Plasma Emission 

Although the first observations of radiation at {ncoc) seem to have been those of Wharton 
[1960], in which noise emission from a pulsed mirror machine was studied, the first extensive 
series of experiments was that of Landauer [1961, 1962, 1964], using a PIG discharge as shown 
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Figure 4. Cyclotron harmonic radiation from a PIG 
discharge [after Landauer, 1961, 1962, 1964]. 



Figure 5. Dispersion characteristics of cyclotron har- 
monic radiation from a plasma [after Tanaka et al., 
1963]. 



in figure 4. These indicated anomalously strong noise components up to the 45th harmonic, and 
fitting the relation {a)/(Oc) = n to better than 2 percent. The most intense radiation was measured 
for the polarization E±B, but harmonic emission was also obtained for E\\B. We may remark 
here that this latter result may easily be predicted from examination of the dielectric tensor obtained 
without making the quasi-static assumption. This indicates that the ordinary wave, E\\B, k±B 
should also show cyclotron harmonic resonances [Stix, 1962; Kubo et al., 1964] . 

The major conclusions of Landauer have been confirmed by Dreicer [1964], in a PIG discharge, 
by Lustig et al. [1964]; Mitani and Kubo [1961]; Tanaka et al. [1963a, b], and Mitani et al. [1964a, b]. 
These latter worked with hot-cathode discharges and were the first to remark that there are signi- 
ficant discrepancies from the relation (a>/a>c) = ^ for maximum noise emission. Typical results 
are given in figure 5. 

The reasons for the very high noise production in PIG discharges are clear from the dis- 
cussion of section 2.5. Such a discharge contains streams of electrons having energy equal to 
the operating voltage, and presents ideal conditions for transverse instabilities or strong beam/ 
plasma interactions with the cyclotron harmonic modes. 

It might be argued that strong directional velocity components were not present in some hot- 
cathode discharge experiments where many harmonics were also observed [Lustig et al., 1964; 
Mitani and Kubo, 1961; Tanaka and Kubo, 1964; Kubo et al., 1964]. It was always noted in these 
studies, however, that noise emission was highest near the cathode end of the tube and at points 
where the plasma entered the magnetic field. The results imply, then, that streams of primary 
electrons entering from the cathode region, and acquiring high transverse velocity components, are 
effective in producing cyclotron harmonic radiation. 

As further support for the hypothesis concerning the effectiveness of charged particle beams 
in exciting cyclotron harmonic resonances, we may refer to radiation measurements carried out 
in the millimeter pressure range where such streams would be rapidly dispersed. Experimental 
results along these lines were reported by Bekefi et al. [1961a, b; Fields et al., 1961; Hirshfield 
and Brown, 1961; Brown and Bekefi, 1962] and have been substantially confirmed by Tanaka and 
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Kubo [1964]. They show characteristic noise radiation dependent on the gas: in hehum and 
hydrogen, the effective radiation temperature corresponds to the electron temperature, but in 
neon, argon, and xenon narrow peaks in the radiation temperature occur, centered on Wc- These 
resuhs can be ascribed to strongly non-Maxwellian electron velocity distributions in the last- 
mentioned gases, where the electron-neutral collision cross section decreases with energy in the 
energy range appropriate to the measurements [Ohara, 1962, 1963], The conclusion to be drawn, 
then, is that beam-type or "humped" velocity distributions are capable of exciting strong cyclotron 
harmonic waves. 

The best test of this, of course, is to use an electron beam to produce the plasma. This ensures 
the presence of both beam and plasma. Two experiments along these lines have been reported 
[Bekefi and Hooper, 1964; Gruber et al., 1964]. In the first of these, the beam and the discharge 
tube were coaxial with the magnetic field, and the gun was only partially immersed, so that about 
15 percent of the longitudinal beam energy was converted into transverse rotation. Emission at 
CDc due to the beam electrons alone could be observed, appropriately Doppler-shifted in frequency. 
When mercury-vapor was present, however, the radiated signal increased by 35 dB and similar 
Doppler-shifted components could be seen up to the fourth harmonic. Radiation from a beam- 
produced plasma in helium up to the 15th harmonic was studied by Gruber et al. [1964], who 
examined the radial and longitudinal amplitude variations of the output signals by means of a 
movable probe. They presented results showing fast-wave growth of the 8th to 10th cyclotron har- 
monics at about 11 dB per wavelength in the direction of the beam, and with rapid radial fall-off. 



3.3. Absorption Studies 

Having discussed experiments in which noise emission only has been studied, we now turn 
to attempts to excite resonance behavior with an external signal. Bekefi et al. [1962] reported 
work on a hot-cathode hehum discharge in a cyhndrical cavity excited in the TEon mode, and im- 
mersed in an axial magnetic field. It was found that for a>^ << w^, resonant absorption occurred 
at o)c, and that for co^ ~ co^, weak absorption of up to the first eight cyclotron harmonics occurred. 

Comparisons of absorption and emission under identical experimental conditions have been 
made by Tanaka and Kubo [1964]. Weak absorption occurred up to the 5th harmonic, whereas 
noise radiation was appreciable up to the 8th harmonic. This is to be expected since we are 
dealing with radiation orders of magnitude higher than that appropriate to thermal equilibrium. 
Reciprocity between emission and absorption would occur only under these special conditions. 

Further absorption studies, both theoretical and experimental, have been made by Buchsbaum 
and Hasegawa [1964] to explain an effect first remarked by Kubo et al. [1964]: the noise radiation 
can show series of subsidiary resonances lying between the main harmonic resonance. With 
increasing discharge current, these successively approach the harmonic resonances, increasing in 
magnitude as they do so, and disappear on reaching them. The subsidiary resonances were ex- 
plained as being due to the presence of radial electron density variation in the plasma. They are 
analogous to the Tonks-Dattner resonances mentioned in section 1, but are "trapped" in the 
center of the column, rather than in the region close to the wall. An alternative source of sub- 
sidiary resonances due to wave interference has been studied recently by Harp [1965]. 

By far the most interesting observations from the point of view of ionospheric propagation 
were the results obtained with the "Alouette" satellite [Calvert and Goe, 1963; Lockwood, 1963]. 
In these, the vehicle carried a transmitter/antenna system giving a pulsed variable frequency sig- 
nal. Strong resonances were found to occur within 1 percent of the calculated local cyclotron 
harmonic frequencies, up to about the 10th cyclotron harmonic. Prolonged ringing was ob- 
served at these frequencies, persisting for many periods after the end of the transmitter pulse. 
Since no such effect has been seen from the ground, we may deduce that the effect was local to 
the transmitter. 

It was suggested that since the perpendicular propagation of cyclotron harmonic modes is 
not Landau damped, and collision frequencies in the ionosphere are relatively low, a ringing 
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effect is to be expected near nojcj where the group velocity is low and the energy is not transported 
rapidly away from the transmitter [Crawford, Kino, and Weiss, 1964a]. This explanation has been 
developed by Fejer and Calvert [1964], and Sturrock [1965], who have estimated the damping 
rate of the resonances. 

3.4. Propagation Studies 

Although the results described so far provide strong support for the cyclotron harmonic wave 
picture, in that they are all consistent with waves propagating nearly normal to the magnetic 
field, one very important observation is lacking. This is a direct measurement of transmission 
across the magnetic field between two antennas. Such an experiment was carried out for the 
first time in 1964 [Crawford, Kino, and Weiss, 1964a, 1964b]. In the original experiment two 
Langmuir probes, parallel to each other and to a static magnetic field, and situated a few millimeters 
apart in a hot-cathode mercury-vapor discharge were used as emitting and receiving antennas 
respectively for a 400 MHz signal. Peaks in the received signal were observed for certain combi- 
nations of discharge current and magnetic field. Typical data are shown in figure 6a. 

A direct comparison with the theoretical curves of figure 3 cannot be made since (Oc and (Op 
(proportional to PJ^) were the experimental variables. It is instructive, however, to compute 
from the theoretical dispersion relation the variation of (co^/co^) with (w^/co^). This is indicated 
in figures 7a, b, c. The coupling of the cyclotron harmonic frequencies and the cold plasma 
upper hybrid resonance are clearly indicated as the parameter (fp/fr) is varied. There is striking 
agreement between the general shapes of these curves and the experimental curves of figures 
5 and 6a. 

If the radial wavelength were constant in the experiment the experimental points could be 
compared directly to the theory. However, the nonuniform density in the cylindrical plasma col- 
umn employed causes the modes to be trapped internally so that the radial wavelength varies with 
current. This is clearly indicated by the presence of two distinct sets of resonances in figure 6a. 
These are for two different trapped modes. When a more uniform discharge is produced by rf 
the transmission resonances are much cleaner. A typical record is shown in figure 6b and indi- 
cates the presence of only one set of peaks. We can, however, indicate on the plot of figure 6a 
the approximate location of the upper hybrid resonance frequency obtained by an independent 
method [Crawford et al., 1963]. This is shown dotted in figure 6a, and reinforces our belief in 
the validity of the cyclotron harmonic wave explanation of the transmission resonances. 
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Figure 6. Experimental observations of perpendicular 
propagation. 

(a) Dispersion characteristics [after Crawford, Kino, and 
Weiss 1964]. 

(b) Transmission record for an rf discharge (the abscissa may 
be calibrated in terms of (coc/o))). 
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Figure 7. Alternative representation of the perpendicular 
propagation dispersion characteristics of figure 3 for a 
Maxwellian transverse velocity distribution [(Vp = (a>/kx), 

VT = (KTe/m)^/2]. 



4. Discussion 



The complicated single-particle analysis required to demonstrate the existence of cyclotron 
harmonic waves has probably been partly responsible for the fact that their existence has been 
established experimentally almost by accident, rather than as the result of deliberate searches. 
Of the studies described, only those concerned with transmission were made specifically to com- 
pare with the form of the theoretical dispersion relation. 

Now that their existence has been established, the importance of these waves can be ex- 
pected to grow. First, there is the possibility of applying them in practical laboratory microwave 
amplifiers and oscillators [Crawford and Kino, 1964]. As they are slow waves, interaction with 
an electron beam directed either along the magnetic field, or across it, will cause spatially growing 
waves. Alternatively, a beam with internal motions can interact directly with slow or fast wave 
circuits to produce gain. Such a beam could be produced in nature by a charged particle stream 
entering a region of increasing magnetic field. Second, taking the opposite point of view, the expla- 
nation of many plasma noise phenomena, including extraterrestrial cyclotron emission, may be 
better understood in terms of cyclotron harmonic wave instabilities [Twiss, 1958; Schneider, 
1959; Bekefi, Hirshfield, and Brown 1961a and 1961b]. A third region of application is in diagnos- 
tics. The "Alouette" ringing experiments suggests that the resonances could be used for meas- 
uring magnetic field accurately. Attenuation measurements could be used to determine collision 
frequencies, and in principle it should be possible to determine electron density and temperature 
from transmission measurements between two antennas. 



801 



All that has been said in this paper about electron effects can be extended to the positive ions. 
Consequently, it is not surprising that experimental observations have already been made of ion 
cyclotron harmonic emission from thermonuclear fusion study plasmas [Bazhanova et al., 1962; 
Pistunovich and Shafranov, 1961; Crawford, 1964c]. The possibility that ion cyclotron waves 
might be useful for plasma heating immediately suggests itself. 



The author wishes to acknowledge the assistance of J. Tataronis with the computations, and \ 

many fruitful discussions of the theory with Dr. G. S. Kino. 

5. Appendix. Integral Representation of the Dispersion Relation for 
Perpendicular Propagation {k^ = 0) 

A convenient integral form for computation may be derived from (17), (20), and (21) by carrying 
out the velocity integrations first. This leads to 

1 = -| exp i^e sin 6 Jo ( 2/Lt sin -j dO (30) 

where d has been written for (coct), and il for (co/coc). Alternatively, this can be expressed as the 
summation 

1 ==-T y exp iQ.e sin 6 Jo (2fjL sin - dO, (Sl) 

Substitution of [(p — 27T{n -\- 1)J for 6 reduces the summation to 

1 =-^ y (exp [- i27ra])«+i exp (ill<^) sin <p A 2fjL sin - d<p. (32) 

If the presence of a small loss component is assumed, this can be further simplified to 

0)2 r2^ ip 

exp [i27Til] — 1 = -f exp (ifiip) sin (f Jo(2/>t sin -^)d(p, (33) 

co^ jo ^ ^ ^ 

which separates into the two forms, 

0)2 r2^ (p 

COS 277(1 — 1 + — I cos Clip sin (p Jo {2^ sin -)d(p 
w^ Jo ^ 

0)2 r 2^ . . ^ 

sin 27rfl = —f sin Clip sin ip Jo (2/x sin -)dip^ (34) 

Wc Jo ^ 

either of which can be used for computation. 

These can be regarded as general formulas from which integral expressions for velocity dis- 
tributions other than the 8-function can be deduced. As an example, for the MaxweUian of (24) 
we must evaluate the weighting integral 

(2^) /o ^" (^'^ '•" 2) "''P (- Sr) 2^''«^"'" (35) 
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The Bessel function identity required is [Bateman Manuscript Project, 1953] 

j ^Jo{2iJi sin |) exp (- Q fidim ^ \ exp [- 2\ sin^ | j . (36) 

The dispersion relation is consequently given by either of the expressions 

cos 27tCL = 1+ -^ exp (— k) cos Cl(p sin (f exp (X cos ip)dip 

0)2 Jo 

0)2 r2^ 

sin27rn=--f exp(— A.) sin flc^ sin </? exp (X cos (p)6/</? • (37) 

o)^ Jo 

The integral relations can be used to determine the series forms of (28) if appropriate Bessel 
function identities are used. For example, differentiation of (18) with respect to /x yields 

sin cp Jo {2^1 sin |) ^ Y — — [/'(aO] sin nS, (38) 

which can be substituted in (34) and integrated term by term. Similarly for the MaxweUian dis- 
tribution we have [Bateman Manuscript Project, 1953] 

00 

exp (X cos (/?) ^ 2 cos mp /,,(X), (39) 

«=o 

which can be differentiated and substituted in (37). Integration, term by term, then gives (28). 
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